Topological entanglement entropy in the second Landau level 
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The entanglement entropy of the incompressible states of a realistic quantum Hall system in the 
second Landau level are studied by direct diagonalization. The subdominant term to the area law, 
the topological entanglement entropy, which is believed to carry information about topologic order 
in the ground state, was extracted for filling factors v — 12/5 and v = 7/3. While it is difficult 
to make strong conclusions about v = 12/5, the v = 7/3 state appears to be very consistent with 
the topological entanglement entropy for the k — 4 Read-Rezayi state. The effect of finite thickness 
corrections to the Coulomb potential used in the direct diagonalization are also systematically 
studied. 

PACS numbers: 03.67.Mn,73.43.Cd, 71.10.Pm 



I. INTRODUCTION 

The physical pictures of states in the first three Lan- 
dau levels of quantum Hall systems are quite distinct. In 
the lowest, n = harmonic oscillator state, Landau level, 
incompressible states are of most importance, typified by 
the Laughlin 1/3 stated and its descendants from either 
the hierarchy^ or the composite fermion^ approaches. On 
the other hand, in the third Landau level (n = 2), com- 
pressible states, stripes, bubbles or possibly liquid crystal 
states are predominate 4 . In many ways, the least under- 
stood behavior occurs in the second (n = 1) Landau level. 
Here, for example, there is an incompressible state at an 
even denominator, v = 5/2^, an incompressible state at 
v = 7/3 reminiscent of the v = 1/3 Laughlin state and 
an odd denominator state, v = 12/5 7,8 , of mysterious 
origin. The present paper is an attempt to better un- 
derstand these filling factors in the second Landau level 
using direct diagonalization and the topological entan- 
glement entropy. In a previous paper—, the topological 
entanglement entropy (the logarithm of the quantum di- 
mension) was used to argue that the essential physics of 
the state at v — 5/2 is in fact given by the Moore-Read 
state. Here a similar analysis is extended to filling factors 
v = 12/5 and v = 7/3. 

The basic idea is that the topological entanglement 
entropy, which of course, is quite distinct from the ordi- 
nary thermodynamic entropy, is an invariant that char- 
acterizes incompressible state o 10 ! 11 . (As a invariant that 
characterizes a certain "isomorphic" class of objects, the 
topological entanglement entropy bears a certain resem- 
blance to the Kolmogorov-Sinai entropy of dynamical 
systems theory. — ) Hence, if the value of the invariant 
for a numerically determined realistic wave function dis- 
agrees with a model wave function, that precludes the 
model wave function. This approach is complementary 
to the calculation of the overlap using direct diagonal- 
ization, in that agreement with a "correct" model wave 
function (or field theory approach) , in principle, improves 
with system size for the topological entanglement en- 
tropy. However, a direct calculation of the overlap gets 



worse for larger systems, even for the "correct" model 
wave function 1 ^. 

The paper is organized as follows: in the following sec- 
tion a few numerical issues are discussed, in particular, 
the role of finite thickness and how to deal with ground 
state degeneracy. In the third section, the results for 
v = 12/5 and v = 7/3 are presented, while conclusions 
are given in the final section. 



II. NUMERICAL CONSIDERATIONS 

The numerical method, developed in ref. [14], [15], 
and notation used in ref. [ 9] are used in the present 
study. The numerical model consists of finite square clus- 
ters with periodic boundary conditions with the electrons 
interacting with long range Coulomb interaction unless 
otherwise noted. The electrons are assumed to be fully 
spin polarized. Direct diagonalization, implemented by 
the Lanczos algorithm 1 ^, is used to find the ground state 
of the system. A brief summary of the numerical proce- 
dure to find the topological entanglement entropy is as 
follows (for full details see ref. [9]). Firstly, for a given 
cluster of N orbitals, the /-orbital entanglement entropies 
are calculated where I is the number of orbitals in the 
subsystem. The Z-orbital entanglement entropy is then 
linearly extrapolated in 1/N to large system sizes; this 
will be termed, "finite size scaling" (FSS). The extrap- 
olated /-orbital entanglement entropies are plotted as a 
function of yl and fit by a line with the y-intercept being 
the negative of the topological entanglement entropy. 

In finding the ground state, the conserved quantities 
are the total number of electrons and the momentum in 
the y-direction. Since we are not doing a full symme- 
try analysis, the observables of the ground state wave 
function are not translationally invariant, even though 
periodic boundary conditions are used. This leads to a 
dependence of the /-orbital entanglement entropy on the 
location of the subsystem. Among the degenerate states 
with different fc y 's a state with equal "left" and "right" 
/-orbital entanglement entropies (see ref. [ 9]) was chosen 
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for which the entanglement entropy was calculated. This 
prescription also works, yields a unique state, at v = 12/5 
and v = 7/3 and this is the choice we make in section 3. 
However, in this section, a different choice will be inves- 
tigated for v = 5/2, for which the previous "recipe" does 
not always yield a unique state. The choice made here to 
guarantee translational invariance of the /-orbital entan- 
glement entropy, for a given cluster size, is to merely av- 
erage the Z-orbital entanglement entropies for the ground 
states with different k y 's. Figure 1 is a graph of the ex- 
trapolated entanglement entropy vs. \fl obtained using 
this procedure. The diamonds were obtained by using the 
pure Coulomb potential and the finite size extrapolation 
(FSS) was done using systems sizes of N e — 12, 13, 14, 15 
and 16 electrons corresponding to N — 24, 26, 28, 30 and 
32 orbitals. Note that the electrons in the n=0 Landau 
level(a spin up and spin down level corresponding to the 
2 in filling factor v = 5/2 = 2 + 1/2) are treated as inert; 
only the electrons in the n — 1 Landau level are consid- 
ered. The topological entanglement entropy, 7, obtained 
from the linear fit was 7 = 2.07 ±0.15. The standard de- 
viation in the individual extrapolated points, indicated 
by the error bars, was always less than 0.14. This is to 
be compared to our previous results, where no averaging 
was done, of 7 = 2.01 ± 0.19 and the maximum error in 
the extrapolated points being 0.18. By either method, 
there is very good agreement with the topological entan- 
glement entropy of the Moore-Read state , i.e. 7 = 2.08. 
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FIG. 1: Entanglement entropy of v = 5/2 state. Entan- 
glement entropy is plotted versus \fl (where I is number of 
orbitals comprising the sub system) for different numbers 
of orbitals. Squares depict results for the system size with 
N e — 16. Crosses depict the same computation, except with 
the more realistic finite thickness Coulomb potential. Finite 
size scaling (FSS), N — > 00, results for entanglement entropy 
versus y/l (diamonds) is shown with linear least squares fit 
(LSF), yielding an intercept of —7 = —2.07 ± 0.15. Circles 
depict latter computation but using finite thickness Coulomb 
potential. 

In figure 1, the squares are the /-orbital entanglement 
entropies for the 16 electron system. Clearly, some sort of 



system size extrapolation is necessary to obtain a value of 
the topological entanglement entropy close to the Moore- 
Read value. An additional feature in the above data (the 
diamonds) is that apart from a linear increase with \fl 
(the area law) there is a superimposed small sinusoidal 
like variation. This is also evident in our previous work 
on v = 1/3 (and other numerical calculations^). We sug- 
gest that this is not an error in the I/TV extrapolation; 
the asymptotic results of ref . [10] and [11] do not pre- 
clude an oscillation which does not grow with the number 
of orbitals. 

Another possibly more realistic choice of interaction, 
which approximately takes the finite thickness of the 
physical system into account was investigated in ref. [18]. 
The Fourier transform of this " square well" potential is 
given by 
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where d is the width of the quantum well. The long 
range (small k) part of this potential is the same as the 
Coulomb potential, the short range (large k) part has 
been modified (softened) by the finite thickness. The 
circles in figure 1 are the extrapolated /-orbital entangle- 
ment entropies and the crosses are the Z-orbital entangle- 
ment entropies for the 16 electron system both calculated 
using this potential with d = 41b, Ib being the magnetic 
length. This choice of d is close to that which gives the 
maximum overlap with the Moore-Read state for small fi- 
nite cluster calculations^. From the (FSS) extrapolated 
entropies, one obtains an topologic entanglement entropy 
of 7 = 2.01 ± 0.19. As indicated by the error bars, the 
extrapolated values have somewhat larger errors (at most 
0.24) for the square well potential. 

As a further comparison of the square well potential 
to the pure Coulomb potential, the previous calculation 
has been extended to N e = 17 electrons in N — 34 or- 
bitals, again averaging over all ground states. These cal- 
culations involve finding the ground states of matrices 
of dimension of order 69 x 10 6 . The results are shown 
in figure 2. The topological entanglement entropy for 
pure Coulomb is 7 = 1.93 ± 0.16 while the square well 
gives 7 = 1.95 ± 0.21. For either choice of potential, 
the results agree less closely with the theoretical value 
of topological entropy for the Moore-Read state. This is 
understandable since the Moore-Read state involves pair- 
ing of composite fermions which is not as favorable when 
the number of electrons is odd and still relatively small. 
The above calculations indicate that as far as the topo- 
logical entanglement entropy and even the entanglement 
entropy is concerned, there is not much difference be- 
tween the square well and ordinary Coulomb interaction. 
Assuming that both potentials give the same phase it is 
not surprising that the topological entanglement entropy 
is the same for both potentials. One could hope that 
the finite size scaling (i.e. the pre asymptotic) behavior 
of one potential or the other would be better. However, 
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FIG. 2: Entanglement entropy of v = 5/2 state including 
N e = 17 data. Crosses and squares are entropy vs. \Tl data 
for N e — 17 with and without finite width Coulomb potential, 
respectively. Circles and diamonds are the N e = 12 — 17 FSS 
extrapolated entropies vs. yfl data, with and without finite 
width Coulomb potential, respectively 

the square well potential does not seem appreciably bet- 
ter and hence forth we have only considered the pure 
Coulomb potential. 

III. TOPOLOGICAL ENTANGLEMENT 
ENTROPY OF v — 12/5 AND v = 7/3 

In this section, the topological entanglement entropy 
of the states at v = 12/5 and v = 7/3 are investigated. 
Due to the possible importance to quantum computing 19 , 
we first consider v = 12/5. One immediately realizes a 
problem in that the largest system size that can be eas- 
ily treated is 14 electrons in 35 states, involving finding 
the ground states of matrices of order 66 x 10 6 . Thus 
only 4 system sizes containing 8,10,12, and 14 electrons 
are included in the extrapolations. The extrapolated l- 
entanglement entropies, as well as the /-entanglement en- 
tropies for 14 electrons are plotted in figure 3. The topo- 
logical entanglement entropy, obtained from a linear fit, 
is 7 = 2.43 ± 0.19. There are particularly large errors 
(standard deviations) in the extrapolated values them- 
selves for I > 8. However, this result seems fairly robust 
in that if one uses only I's up to 8 to get the linear fit, 
one still gets a topological entanglement entropy of 2.43 
with a larger standard deviation of 0.31. Taking into ac- 
count the two boundaries for our system on the torus, the 
k = 3 Read-Rezayi state has a topological entanglement 
entropy of 2 In .D « 2.90 where D = 2 . „ . Hence, for 
our value of 7 » 2.43 to be consistent with k = 3 Read- 
Rezayi one needs to take the large possible errors in the 
extrapolated values into account. Recall that this was 
not necessary at v = 5/2 with the Moore-Read state and 
the extrapolated errors in the /-entanglement entropies 
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FIG. 3: Entanglement entropy of v = 12/5 state. Diamond 
symbols depict entropy vs. y/l data for N e = 14 electrons in 
N = 35 orbitals. FSS extrapolation of entropies from N e = 
8— 14 yields entropy vs. \fl data, depicted with circles. Linear 
least squares fit (solid line) gives an intercept of —7 = — 2.43± 
0.19 for the v = 12/5 state. 



were smaller. 

From an experimental standpoint, the state at v = 7/3 
is more robust then that at v = 12/5 7 . Hence there is 
some hope that numerically it will be easier to study 
the state at v = 7/3. In figure 4, we plot, as the cir- 
cles the extrapolated /-entanglement entropy vs. y/l for 
an extrapolation involving 9,10,11 and 12 electrons; the 
diamonds are an extrapolation for 10,11 and 12 elec- 
trons. The squares are the 1-entanglement entropies for 
the largest system size, 12 electrons. A linear fit to the 9- 
12 electron extrapolations gives an entanglement entropy 
of 2.92 ±0.24, with rather large error bars on the last four 
extrapolated /-entanglement entropies. By way of com- 
parison, the topological entanglement entropy expected 
from the Laughlin state at v — 1/3 is approximately 1.09 
and for the k = 4 Read-Rezayi state is 21nZ? « 3.58 
where D = » 6 m . Since the error bars in the 10-12 ex- 

2 sin 

trapolations are smaller, it is perhaps more meaningful 
to do the linear fit in this case; one finds a topological 
entanglement entropy of 7 = 3.56 ± 0.33. This is in very 
good agreement with the value expected for the k = 4 
Read-Rezayi state. 



IV. CONCLUSION 

We have calculated the topological entanglement en- 
tropies at v = 12/5 and v = 7/3 in the flat torus ge- 
ometry. There are large possible errors associated with 
our results, due to the difficulties in extrapolating to 
large system sizes. That being said, still our calcula- 
tions tend to support the picture that v = 12/5 corre- 
sponds to the k — 3 Read-Rezayi state(numerical evi- 
dence for this is presented in ref. [20]) and v = 7/3 
corresponds to k = 4 Read-Rezayi (numerical evidence 
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FIG. 4: Entanglement entropy of v = 7/3 state. Same as 
figure 3 except: squares are N e = 12; diamonds represent 
FSS extrapolation for N e — 10 — 12; circles represent FSS 
extrapolation for iV e = 9 — 12. Linear least squares fit is 
performed for N e = 10 — 12 data (solid line) yielding 7 = 
3.56 ± 0.33 for the v = 7/3 state. 

that v — 7/3 is not a Laughlin state is given in ref. [21]) . 
Qualitatively, we find that the topological entanglement 
entropy is larger for v = 7/3 then for v = 12/5, this 
is consistent with assigning these fillings to the k = 4 
and k = 3 Read -Rezayi states. The quantitative case 
for the Read- Rezayi states is stronger for v = 7/3 in 
that one can get agreement without taking the error 
bars of the extrapolation into account. Taken at face 
value, there is excellent agreement between our numer- 
ical result for v = 7/3, topological entanglement en- 
tropy of 3.56 ± 0.31 and the value for the k = 4 Read- 
Rezayi state « 3.58. Other theoretical alternatives, the 



Laughlin state or the hierarchical construction (based on 
v = 5/2) of Bonderson and Slingerson™ (topological en- 
tanglement entropy of 2 In (2\/3) ~ 2.482^) have lower 
topological entanglement entropies which are not consis- 
tent with our numerical results. The v — 12/5 of Bon- 
derson and Slingerland has topological entanglement en- 
tropy of 2^(2^7(5)) « 3.0Q22. This is consistent with 
our numerical results ( or rather as consistent as the k=3 
Read- Rezayi state). Qualitatively however, this disagrees 
with the numerics in that the topological entanglement 
entropies of the Bonderson-Slingerson states increase in 
going from v = 7/3 to v = 12/5. 

If the k — 4 Read-Rezayi state captures the physics of 
v = 7/3 we are left with the following paradoxical situa- 
tion. Experiment indicates that the particle-hole conju- 
gate of v = 7/3, i.e. v = 8/3 has a quasi particle charge 
of e* = e/324. Of course, if particle- hole symmetry is 
not badly broken this leads to a quasiparticle charge at 
v = 7/3 of e* = e/3. However, the quasiparticle charge 
of k = 4 Read-Rezayi is e/6. Due to the possibility of 
breaking particle-hole symmetry there is not necessarily 
a conflict with experiment and the lack of a incompress- 
ible state at v = 13/5 7 indicates such symmetry breaking 
is present. In any case, a direct experimental measure- 
ment of the quasiparticle charge at v = 7/3 would be 
very interesting. 
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